Conditions for plasma evolution to the strong general 

Woltjer- Taylor state * 

Jun-hua Chen 1 ^and Hong-yi Fan 1 

1 Department of Material Science and Engineering, 
University of Science and Technology of China, 
Hefei, Anhui, 230026, China 

December 27, 2012 



Abstract 

We find that the proof in the recent paper ^ can not justify the authors' conclusion. We provide a real 
proof that any state will eventually evolves to the Woltjer- Taylor state exponentially. However, this kind 
of evolution is is mainly due to Joule heat, which also makes the magnetic field vanishes exponentially. 
Zero Woltjer- Taylor states are not physically attractive. Instead of examine A, we introduce the quantity 
9 Vx j j and R to examine if the plasma reaches to the strong (general) Woltjer- Taylor state, and then 
derive the condition for the evolution to the strong/general Woltjer- Taylor state. 
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1 Introduction 



H = J v A-BBd 3 x, 



Q = J v A 2 d 3 x, and B is the mag- 



Based on the strong astrophysical and laboratory 
evidence that plasmas tend to evolve towards the 
Woltjer- Taylor state (WTS) Etini, characterized by 
V x B — aB, and being unsatisfied with Tay- only if A = 
lor's theory (conjecture) EHI3I ) Q m et al ^ de- 
veloped a theory which they claimed in favor of 
the evolution towards to the WTS. They introduced 
A = QW — H 2 ^ for measuring the deviation 



netic field, A is the vector potential. They pointed 
out that equality A = holds if and only if B = aA 
for some constant a and therefore VxB = aB. Qin 



et al then proved that ^A ^ and gjA = if and 
Based on the above inferences they 
drew their conclusion that lim A (t) =0 and the 

t— >+oo 

WTS is reached. 



of the plasma to the WTS, where W — J y B 2 d 3 Q 



However, the following counterexample shows that 
Qin's conclusion is not justified by their argument if 
there is no further improvement on the estimation of 



*Work supported by the National Natural Science 
Foundation of China under grant 11105133 and 11175113 

i~To whom correspondence should be addressed. Email: 
cjh@ustc.edu.cn 



Let x (i) ^ be a function satisfying ^ 



We see that < 0, and 



di 



if and only 



if x — 0. The solution to the equation % 



T+F 



1 



x(t) = x (0) exp (- arctani). Though ^ < for all 
t, we still have x (+00) = x (0) exp (— §) > 0. 

Further, even if one do manage to prove that 

lim A (t) — 0, since A is quartic in B, A's ap- 

t— ¥+00 

proaching zero could possibly be the simple result of 
the vanishing of B, as we will see later in this paper. 
Instead of examining A as Qin et al did, we intro- 
duce a dimensionless quantity, i.e., the angle Q^ x § § 

between two fields B and VxB, which is defined by 



cos 9 



V x B,B 



VxB.B 



V x 5, V x B) (B,B 



(1) 



As one can see, Vx g g measures the deviation of 
plasma to the strong WTS ( by "strong" we mean 
the case of V x B = aB and B ^ 0, a is a constant 
). Further, if a (x) in equation V x B = aB varies in 
space, we say the plasma is in general WTS. In this 
case we can introduce another dimensionless quan- 
tity, the relative residual error R (see Eq. (|2T]) ). As 
we shall explain below Eq. ([2"T]) . R measures the de- 
viation of plasma to the strong general WTS. By cal- 
culating the evolution of 0y x g § and R, we can check 
whether plasma evolves to strong (general) WTS. 



2 Presetting 

For physical simplicity, we will use SI unit and 
Coulomb gauge V • A — in this paper. The scalar 

3^?/ 



Pc{x ,t) 



'■^-Lrrdrx 



potential takes the form ip(x,t) = J g^ng^gri ' 
under Coulomb gauge condition. Since there is no 
net free charge p e in plasma, we have if (x, t) = 0. 

The vector potential A (x, t) related to a uniform 
magnetic field Bq (t) is unbounded in space, therefore 
we can not expand such A into Fourier series. Fortu- 
nately, uniform magnetic field Bq (in full space, or in 
a sufficiently large volume V) does not exist in real 
physical world. Moreover, we can safely assume that 
A and B both vanish at the boundary dV of V. In 
the following discussion, we assume that all the fields 
A, B, E, J, v are "good enough" so that they can be 
expanded into Fourier series as 



A, B, E, J, v) = £ (4, 4, Jg, 

k^0 



'ik-x 



(2) 



For a finite volume V, the magnitude of nonzero 
k in Eq. ^ has a universal lower limit ko > 0. fco 
depends on the size and shape of V and the boundary 
conditions only. For example, let V be a cubic with 
side length L and choose the boundary conditions so 
that the fields vanishes at the boundary, then fco = j- ■ 

The inner product (^A, Yj of two real fields X, Y 

in V is defined as (^A, Y^ = J v X-Yd 3 x. For any real 

field A wc have ^A^ = X The inner product 
becomes 



\fc#o / \fyo J 

= VY J X r Yi = VY J X* r Y % . 



(3) 



By Cauchy-Schwartz inequality, y = 

(A,A)(f,f) - (X,Y) 2 > 0, = if 



and only if A = aY , where a is a constant. There- 
fore &x Y does measure the deviation of the two 
fields to the equation A = aY. 



3 Evolution of the Magnetic 
Field 

We have the following equations in real space and in 
the Fourier component space 



V-E = 0, k-E % = Q. 



V-B = 0, k-B r = Q. 



Taking time derivative on Eq. ([5]), we have 



JTl 

k-B^ =k-—B r = 0. 
k dt n k 



(4) 
(5) 

(0) 
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Equations B = V x A and V • A = give 

ik x B r 



k 2 



(7) 



Faraday's law reads 



OB d 
V x E = — —,ik x Er = — -Br. 

at k dt k 



By Eq. 0, we can solve from Eq. © 

iife x 

p at k 

^k — fc2 

We also have Ohm's law 

J = a{E + vxB], 



(9) 



(10) 



where a is the conductivity of the plasma. For sim- 
plicity we assume that a is a constant, as what Qin 
et al did ^ . Let A and r be the typical scale of field 
variation in space and time respectively, in the case 
that v <C ^, we have 



iak x 4: Br 
~ k 2 



Equation V x B = fi J + e Mo^f gives 
ik x Br ikx £>B*r. 



Mo 



'jfe i _ v,v ^ dt 2 k 

£Q k 2 



(11) 



(12) 



When the fields are varying slowly, i.e., ^ < c, we 
have 

ik x Br 

J %™ — " — • (13) 
Mo 



Combining Eq. ([6]). (fTT]) and (JTSJ) we have 

d ~ k 2 
-Br « Sr. 



(14) 



Therefore 



^(t) = B jg (0)exp(-— t 



(15) 



We see that Br (t) goes to zero exponentially. 
Br (i) with larger k 2 varnishes faster. When t — > +oo, 
(8) Br with minimum k dominates. The fine spatial 
wrinkles (short wavelength fluctuations) of B will be 
ironed out by Joule heat exponentially. 



4 Conditions for Approaching 
Strong WTS 

Now we are ready to calculate the evolution of 
9 vxbJ- Firstly, 



V x B,V x B) = V J2 

k^Q 

= vj2k 2 s fe -(o) 12 



ik x Br 



2k 2 - 



exp I —-=—t 



(16) 



fc#0 



2 / 2k 2 

exp 1 ) , (17) 

V ^Mo 



V x B,B) = V J2 ikxBr(0) • 55 (Q) e~^5\ 

k^O 

(18) 

Let fc m ; n ^ fco > be the minimum value of k 
with Br (0) 0. When i — ► +oo, terms with A; 

fc min dominate in Eqs. (UHl HZ1 HE]) ■ We have 



lim cos I 

t— >+oo 



VxB,B 



(19) 



E ikxB % (0) -BS(0) 

|fc|=fcmin 



|fc|=fcmin 



E %(o) 



COS (5, 
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which is a constant depending on the initial spa- 
tial distribution of B only. If sin S ^ 0, plasma 
can never evolve to strong WTS. More precisely, we 

2/c j. 2 /c 



B = in V, which means V x B = a (x) B. 
We have 



have cos / 



VxB.B 



cos S + O 



exp 



where k nex t is the minimum value of 



> &min with 



the rate O 



exp i 



afia 



if sin 5 = 0. Re- 



(V x Bj 

therefore 
E 



x B 



E rx 4-) x%-p ( 22 ) 



calling the fact that Bj (t) = B £ (0) exp ( -^*J < 

B^ (0) exp ^— ^J^^j j we see that a reasonable re- 
quirement is k^ ext > 2k^ ain so that plasma can ap- 
proach strong WTS "before" B practically goes to 
zero in the case sin S = 0. 

As a by-product of Eq. (fl~5)). we calculate 

A = ( A. A) ( B. U) - ( . 1 . S ' 2 



i? 



2 _ fc ^° 



S (<fx5 r (0))xBg_ r (0)e ""0 



E 



E B E (0)-BI(0)e -"0 



E (fxB r (0))xBg_ r (0) 

|f|=z|fc-f|=fc min 



E 

.MO 



fc 2 



E 



E %(o) 



%(0) 



which does decreases exponentially with time at the - 2 jy V • (b x eJ d a x 2 (e. VxB 

exp ^ ■ But this is merely a triv- = 2 J av , (^B x E^j ■ dS — 2 ^B, /io J + eo/io 



fci.. I E -B £ (o)-b:(o) 

v l«l=*mln 

(23) 

as t — ¥ oc, which is also a constant and nonzero in 
general cases. 

5 Further Discussions 

In the most general cases, W will decay at the rate 



f =2 f,B) = -2(Vx£,B 



dt 



di 



rate O 



ial inference of t 



O 



exp 



ie fact that B vanishes at the rate 



If a in equation V x B = oB is not a constant, as 
is widely discussed, the plasma is in general WTS. In 
this case sin6* Vxi j g is generally nonzero, therefore 
®VxB S can n °t measure the deviation of plasma to 
the general WTS. A better quantity that measures 
the deviation of plasma to the general WTS is the 
relative residual error R defined by 



f)E 
dt 

-2e 0f x -2hq(e J 

2/i (uxB,/>-^J,J 

cr max 

J24) 

where (T max = sup [a (x)] < oo. We see that B will 



varnish due to Joule heat within time scale 



for 



R 2 = 



VxB x B, V x B xB 



(21) 



It is obvious that R = if and only if ( V x B ) x 



any real plasma with finite conductivity a (x) . The 
adjective "strong" is necessary when talking about 
WTS in physics, and the corresponding time scale is 
indispensable when talking about the evolution to- 
wards WTS, since lim B = 0. 

t— >+oc 

In the above discussions we see that under the as- 
sumptions we have made, the WTS is reached only if 



4 



the initial distribution of B satisfies some constrains. 
But astrophysical and laboratory observations indi- 
cate that WTS is the general status for plasma. The 
only significant factor that we have not taken into 
account above is the non-homogeneousity of a. 
For real plasma with inhomogeneous a — oq + 



E ^ lk - s , 

k^O 

i(k-T) 



Eq. (fTl) becomes Jt ~ — Co 



ikx4tBt. 



(*-0 



Therefore Eq. (fT4|) becomes 



— R 

dt D k 



E 



Tr kx 



k z B t 
CT oMo 



(25) 



(*-9 



Plug the first order approximation Br, (t) » 
Bj: (0) e "ofo in the last term in Eq. (|25|) . we have 



dt k &Qfl>Q 

<Jr fex[(fe-f)xB £ _ r (0)] _(£dl! t 

\ _! Li I £_i ig o- 0( i 



(26) 



We see that even for k > /c m i n , there could still ex- 



ist some ( ^ so that 



k-l 



will generate one term proportional to e '«™ 



dominates in Br. when t —> +oo and change the val 



= &min < k, which 

that 
val- 

Thus the non- 



ues of lim L v s s and lim R 

homogeneousity of a is one very possible candidate 
for the real power that pushes plasma towards (gen- 
eral) WTS before B practically goes to zero. 
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